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Abstract. We present in this paper a routine which construct the ideal generated by a list 
of elements in a matrix Lie algebra at any particular characteristic. We have used this 
algorithm to analyze the problem of the simplicity of some Lie algebras (see the refer- 
ences). Notice that it is not included the particular construction of a basis (chosen a Lie 
algebra) denoted by {xj }" =1 . 



■ Previous definitions 

As we have comment in the abstract, there are some commands that must be pre-definedin order to use 
the algorithms, listed below: 

1. A basis of the Lie algebra, denoted by x, , ( = 1, n. 

2. The dimension of the Lie algebra, denoted by n . 

3. A function which recognize a element of the Lie algebra, denoted by Recog, i.e. write an element as a 
non-evaluatedlinear combination of the elements of the basis. 

As an example, we present the Lie algebra M 2 (K) with char(7T) = k. Then we define 
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Xl = (o o) ;X2 = (o 5) ;X3 = (? o) ;X4 = (o J)'' n = 4 ' 

Recog [m_] : = Module [ { } , 

Ai =m[[l, 1]]; A 2 =m[[l, 2]]; A 3 =m[[2, 1]]; 
A 4 = m[ [2, 2] ] ; Sum[Ai * (HoldForm [x] ) ± , {i, n}] ] 

For instance, we make 

R6COg [(c d)] 

ax\ +bx2 + c X3 +dx4 

It is also important to remark that we have not included any filter in order to check if the data is well 
created. 

■ Routines 

We define the Lie bracket 

c[x_, y_] : = x.y - y.x 

and we can, for example, construct the multiplication table of the Lie algebra 

Table [Recog [c [xi , Xj ] ] , {i, n} , {j, n}] 
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The function 



NonZero [lis_] := Module [ {new ={}} , 
Do[ 

If [Not [lis [ [i] ] = = = * lis[ [i] ] ] , AppendTo [new, lis [ [i] ] ] ] , 
{i, Length [lis] }] ; new] 

works as 



NonZero[{l, 2, 0, 3, 0, 0, 4, 5, 0}] 

{1,2, 3,4,51 



The following comands 

Coord[y_] : = Module [ {coo, aa} , coo = Array [aa, {n}]; 

Do[aa[i] = Coefficient [y, (HoldForm [x] ) ± ] , {i, n}]; 
coo] 

Comb [ lis_] : = Module [ { } , 

Sum[lis[[i]] * (HoldForm[x]) i , {i, n}]] 

can be explained with the following two examples 
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Coord[Recog[(* J]] 

{a, b, c, d] 



Comb[{l, 2, a, 4}] 

X\ +2x2 + 0x3 +4x4 

We have then the principal algorithm. It needs two entries: a list of element (perhaps only one) and the 
characteristic of the base field. As a result, we have the ideal generated by the list of elements, a basis, 
and a list of partial results in terms of the depth of the Lie brackets. 

IdealGenerated[lis_, char_] :=Module[{k, ideal}, 

ideal = Map [Comb, NonZero[ 
RowReduce [Map [Coord, Map[Recog, lis] ] , Modulus -> char] ] ] ; 
Print["Depth = ", 0, " -> ", ideal]; 

dim[0] = 0; dim[l] = 1; k = 1; 

While [ 

dim[k] < n && dim[k] # dim[k - 1] , conj = ideal; k = k + 1; 
Do [AppendTo [ ideal , Recog [c [ReleaseHold[con j [ [ j] ] ] , x ± ] ] ] , 
{i, n}, {j, Length [conj] }] ; 

ideal = Map [Comb, NonZero[ 
RowReduce [Map [Coord, ideal] , Modulus -> char] ] ] ; 
Print ["Depth = ", k-1, " -> ", ideal]; 
dim[k] = Length [ideal] ] ; 

Print["Ideal <", Map[Recog, lis], 
"> = ", ideal, " with ", " dimension = ", 
dim[k], " and char(K)=", char]] 

Using our example, we can compute 

IdealGenerated[ {x 2 } , 2] 

Depth = -> [x 2 1 

Depth = 1 -> [xi + x 4 , x 2 } 

Depth = 2 -> [x\ + x 4 , x 2 } 

Ideal <{x 2 }> — {X\ + %4 , X2 } with dimension = 2 and char(K)=2 

IdealGenerated[{xi } , 2] 

Depth = -> \x x 1 
Depth = 1 -> \x { , x 2 , x 3 } 
Depth = 2 -> \x { , x 2 , x } , x 4 } 

Ideal <{ j with dimension = 4 and char(K)=2 

IdealGenerated[ {xi , x 2 } , 2] 

Depth = -> {x x , x 2 } 
Depth = 1 -> {*! , x 2 , x 3 , x 4 } 

Ideal <{xj , x 2 }> = {x { , x 2 , x 3 , x 4 } with dimension = 4 and char(K)=2 
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IdealGenerated[ {x 3 , x 3 - xi } , 3] 

Depth = -> , x 3 } 
Depth = 1 -> {x\ , x 2 , x 3 , x 4 ) 

Ideal <{x 3 ,x 3 - X\ }> = [x\ , x 2 ,x 3 , x 4 ) with dimension = 4 and char(K)=3 

IdealGenerated[ {x 2 } , 3] 

Depth = -> [x 2 1 
Depth = 1 ->{x\ + 2 x 4 , x 2 } 
Depth = 2 -> {xi + 2 x 4 , x 2 , x 3 } 
Depth = 3 -> [x\ + 2 x 4 , x 2 , x 3 } 

Ideal <jA: 2 }> — {X\ + 1 X4 , X2 , X3 j with dimension = 3 and char(K)=3 
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